Vibration-induced Kondo tunneling through metal-organic complexes with even 

electron occupation number 
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We investigate transport through a mononuclear transition-metal complex with strong tunnel 
coupling to two electrodes. The ground state of this molecule is a singlet while the first excited 
state is a triplet. We show that a modulation of the tunnel-barrier due to a molecular distortion 
which couples to the tunneling induces a Kondo-effect, provided the discrete vibrational energy 
compensates the singlet/triplet gap. We discuss the single-phonon and two-phonon assisted co- 
tunneling and possible experimental realization of the theory. 

PACS numbers: 85.65.-|-h, 73.23.Hk, 73.63.Kv, 63.22.-|-m 
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The effect of quantized vibrational motion of a single- 
molecule device on the electron transport has been inves- 
tigated intensely. Vibrational effects have been observed 
in the sequential tunneling regime P, Qj 13 well as 
the strong tunneling regime 0,11, . Theoretical atten- 
tion first focused on the weak coupling limit, where ex- 
change of mechanical energy quanta with tunneling elec- 
trons (vibration assisted tunneling) 0, 0, S E 
and modulation of the tunnel barriers (shuttling) 
were discussed. In the strong tunneling limit the effect of 
vibrations on the Kondo-anomaly in the linear conduc- 
tance was discussed due to assisted tunnelin g ll^ lisl 
and also due to tunnel-barrier modulation In this 

Letter, however, we demonstrate that discrete vibra- 
tions through tunnel-barrier modulation can induce a 
Kondo effect in a transition-metal (TM) organic complex 
(TMOC) with even number of electrons by compensat- 
ing for the singlet-triplet splitting at zero bias. Already 
for molecules of moderate size, many vibrational modes 
are available, underlining the importance of considering 
the above effect. Our effect is essent ially different from 
the phonon-assisted Kondo effect in [l^ [13 : no strong 
electron-vibration coupling is required nor special gate- 
voltage restrictions. 

We study the transport through a TMOC with a TM ion 
secluded in a ligand cage. The cage is in tunnel contact 
with metallic reservoirs (surface, STM nanotip, or edges 
of metallic wire in electro-migration or break junction 
geometry). Fig. ^a) illustrates this setup. We consider 
a TMOC with even electron number TV fixed by charge 
and energy quantization. The ground state is supposed 
to be a spin singlet, and the energy of the lowest triplet 
excitation A exceeds Kondo temperature Tr- . The linear 
conductance is thus suppressed. To investigate how in- 
tramolecular vibrations may induce transport through a 
Kondo effect, in the first place one should incorporate a 
vibronic mode in a generic tunneling Hamiltonian 
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FIG. 1: Schematic situation: (a) Electrode tunnel-coupled to 
a transition-metal organic complex. Charging of the complex 
by a tunnel process deforms the outer part of the ligand cage 
without strongly afi'ecting the direct coordination-sphere of 
the metal ion and thereby the ligand-field splitting. We as- 
sume that the extra electron is localized mainly on the cage. 
Electrons tunnel onto the ion through the tails of the molec- 
ular state centered on the ion, which includes admixtures of 
the outer shell electronic states. Therefore the main effect of 
the charging is the modulation of the tunnel barrier between 
the ion-centered states and electrode, (b) Electronic e-type 
states discussed in the text in order of increasing energy. 



Here H,noi includes the 3d electron levels in a ligand 
field of the cage electrons, the molecular orbitals of these 
ligands, as well as interactions within the 3d shell and 
within the cage. One should take into account the three 
most relevant charge states including their dependence 
on the vibrational coordinate of the cage Q: 
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The last term T„ is the kinetic energy of the cage distor- 
tion. The eigenstates of Hq''^^^ are admixed to those of 
iy(^) by the tunneling Htun of electrons from the reser- 
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tunneling Hamiltonian may be rewritten in the form 
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H.c. 



(3) 

Effectively a single electrode remains after the standard 
rotation of electron states jl^l- We assume that there 
is no direct tunneling contact between the 3d electrons 
of the caged TM and the electrodes. However the Zd- 
orbitals are hybridized with ligand molecular orbitals, 
and tunneling becomes possible due to small overlap be- 
tween the "tails" of distorted orbitals d and the lead elec- 
trons (see caption to Fig. ^a) for more details) . 

We write Hmoi for the low-lying many particle states 
of the TMOC in the compact form 
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where X^^' = |A)(A'| are so called Hubbard operators 
describing transitions between the eigenstates A, A' of 
the TMOC, and all many-particle effects and their co- 
ordinate dependence are included in the energies Ea{Q). 
Similarly, states of the charged TMOC with iV -> ± 1 
are denoted as I7). The microscopic origin of the low- 
energy singlet (S) and triplet (T) is most easily conceived 
as follows. The evenly occupied ligand cage is in a sin- 
glet spin state, and it does not affect the structure of 
the spin multiplet. We assume that the ligand field of 
the cage has low symmetry, so that the 5-fold orbital 
degeneracy is completely lifted. For instance, for a dis- 
torted tetrahedral symmetry of the ligand field we fo- 
cus on the configuration d^(e^) with two e-states split 
by S due to this distortion. In case of distorted cubic 
symmetry, the same S/T multiplet arises for the config- 
uration d^{t^e'^). Fig. njb) illustrates the states of the 
TMOC: it is seen that in the case of a weak intrashell ex- 
change I < S the ground state is singlet S and the lowest 
excitation is the S/T transition. The energy difference 
A = Et — Es = (5 — / is assumed to be larger than the 
Kondo temperature Tx- The excitation energies of two 
singlet excited states S' and S are S and 26, respectively. 
We ignore below these singlet excitations, since they are 
not involved in the Kondo tunneling. 

We use the simplest approximation of single harmonic 
vibration mode with frequency Q. We assume that the 
relative shifts of the harmonic potentials in different ex- 
cited and charged states is negligibly small i.e. for all 
X = A, 7 we have E^{Q) « E^ + nQ'^/2 and T„ = 
The coordinate Q is normalized to the zero-point mo- 
tion. In this weakly non-adiabatic limit vibrations are 
slow compared to the electron motion on the molecule. 
No vibrational excitations can be induced by the tunnel- 
ing in the linear transport regime, unless the distortion 
of the ligand cage Q affects the tunnel amplitude w{Q). 
Then local phonons can be emitted or absorbed. The 



Htun = W{Q) J2 J2 [^^^^ka + H.c] (4) 
k A-ycr 

(here the second summation is restricted by spin selection 
rules). This situation is opposite to the anti-adiabatic 
limit where the vibrations are fast and many phonon ex- 
citations are involved in the Kondo exchange [l^ . 
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FIG. 2; Two types of phonon-assisted Kondo co-tunneling 
processes, (a) Virtual phonon absorption initiates a S/T tran- 
sition, Kondo processes take place in an intermediate triplet 
state and the phonon is emitted in the end. (b) Every spin 
flip process in the intermediate triplet state is accompanied 
by a two-phonon process. Points and crosses denote spin-flip 
process in the S/T and T/T channels, tij are initial and final 
times, and ti...t„ denote the intermediate co-tunneling acts. 

Thus we consider the effect of the modulation of the 
tunnel amplitude on the linear transport in the Coulomb 
blockade regime. In this regime one has to eliminate 
the tunneling term from the Hamiltonian by summing 
over virtual processes where A^ ± 1 electrons occupy the 
molecule and vibrational quanta are excited in these vir- 
tual states. This procedure known as a Schrieffer- Wolff 
(SW) transformation, leads us to a Kondo Hamiltonian 
(Eq. ® below) for a S/T multiplet [13 and an oscillator 
with electron-phonon interaction built into the exchange 
coupling. According to Ref. [T^ . transitions within the 
singlet/triplet spin manifold are described by two vectors 
S and R, where S is the usual spin 1 vector and R is a 
vector describing S/T transitions. These two vectors are 
constructed by means of Hubbard operators X 
following way [T^ : 



in the 
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The effective Hamiltonian arising after the SW trans- 
formation has the form 



H, 



eff = Hres + ^AS^ + J^S • S + JflR • S + ^ (6) 



The electron spin operator is given by the conventional 
expansion s = ^ Zlfefe' Zlo-o-' ^L'^cr'j'C/c'o-' where x is a 
Pauli vector. The exchange coupling constants Js.b{Q) 
are estimated as Js{Q) ~ J2y — Ej\, and 

Jr = aJs- Here a < 1 is a coefficient arising because of 
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admixture of singlet states 5", S" to the ground state [l^ . 
In the weakly non-adiabatic regime the kinetic energy in 
the denominator has been neglected and the nearly iden- 
tical Q dependence are practically cancelled out in the 
addition energies. Thus, no multiphonon replicas appear 
in the denominators of Js,b. unlike Ref. [14]. 
The main source of phonon emission/ absorption in our 
case is the tunneling rate |w((5)p. Expanding it in the 
quantized displacement operator Q — (b^ + b)/^/2 we 
come to phonon assisted exchange vertices presented in 
Fig. 131 In accordance with Fig. [3 we retain only single- 
phonon processes for Jr{Q) and only two-phonon pro- 
cesses for Js{Q): 

JsiQ) = Js + jsQ^ MQ) = Jr + jrQ- (7) 

It is obvious that js <C jn- To draw these vertices we 
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FIG. 3; (a): Bare exchange vertices Js,r', (b) Single phonon 
correction ju to the vertex Jr; (c) Two-phonon correction js 
to the vertex Js; (d,e) Renormalized vertices 71,2 correspond- 
ing to the processes illustrated by Fig. ^a), (b), respectively. 

used the fermionic representation for the operators lO 
= \/2(/oV-i + ./iVo), = flh - (8) 



V2{flfs - = -(/oV. + /]/o). 



The spin-fermion, electron and phonon propagators are 
presented in Fig. El by solid, dashed and wavy lines, re- 
spectively. Physically, Figs ^a.) and ^d) describe the 
regime, where the virtual phonon absorbed by the TMOC 
in the process of co-tunneling brings the energy neces- 
sary to compensate the S/T gap A. Then the Kondo 
effect develops as multiple spin-flip co-tunneling in the 
triplet state, and the system returns to the ground state 
singlet after virtual phonon emission. The vertex correc- 
tions are calculated by means of analytical continuation 
of Matsubara-type diagrams from imaginary axis to the 
real frequency axis. Summation of all parquet diagrams 
entering the 4-tail box in Fig. Efd) gives for 71 
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here A ~ 1 is a constant determined by spin algebra. The 
Kondo temperature extracted from this equation reads 
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Here D is the effective width of the electron conduction 
band and p is the density of states on the Fermi level. The 
second channel illustrated by Figs|2Ib) and|2Ie) involves 
two single-phonon processes compensating S/T transi- 
tions. The exchange acting in the intermediate triplet 
states is accompanied by too-phonon processes. Then 
summation of parquet diagrams entering the 4-tail ver- 
tex in Fig. ISJe) gives for 72 
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The Kondo temperature characterizing this channel de- 
pends on the phonon-assisted exchange constant js 
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One concludes from these calculations that the single- 
phonon processes are sufficient to compensate the en- 
ergy of the S/T splitting and induce resonance tunneling 
through the TMOC provided a local vibration mode with 
appropriate frequency satisfying the condition 



\n - Al < T; 



(1) 
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(13) 



exists in the cage. One can expect in this case a signifi- 
cant enhancement of the tunnel conductance already at 



T > T^^^ according to the law G/Gq ~ In" 



'■{T/Tj, 



(1) 



where Gq is the conductance at unitarity limit T ^ 0. 

Thus we formulated in this paper conditions under 
which phonons are not only involved in Kondo screening 
but even induce Kondo tunneling. In is worth mention- 
ing that in spite of the fact that the Kondo effect exists in 
our case only under phonon assistance^ the Kondo tem- 
perature (|12|l is the same as in the usual Kondo effect. 
Since Tfj- is high enough (~ 10 K) in electro-migrated 
junction experiments with a TMOC deposited between 
contacts 0,11, 0] , the effect predicted in this work seems 
to be easily observable. The crucial point is the existence 
of phonon satisfying condition 1)13(1 in a TMOC with the 
S/T multiplet as a lowest spin excitation. One should 
note, however, that even if this condition is not exactly 
satisfied, one may tune the system by applying the mag- 
netic field. Then the triplet is split, and only the level 
Ex I = Et — Ez is involved in the phonon induced Kondo 
tunneling [Ez is the Zeeman energy). Then A in 1(1 3|l 
is substituted by A^ = A — Ez, and Ez may be tuned 
until the inequality is satisfied. Thus the vibration gives 
rise to a magnetic field induced Kondo effect at Zeeman 
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energies which can be much smaller that A. The Kondo 
screening takes place due to the processes presented in 
Figs |2Ia) and|3fd). The only difference is that in this 
case the effective spin of the TMOC is one half instead 
of one • The theory of this effect will be presented in 
a separate publication. 

Another way to tune the condition 113() is to stretch 
the break junction and thereby distort the TMOC and 
change the frequency fl. Such mechanical control was 
demonstrated experimentally for the H2 molecule [25j . 
There it was also shown that the isotope effect may be 
used for the same purpose. 

The generic feature of the phonon-assisted Kondo 
screening discussed in this work is that only the virtual 
phonon excitations enter the co-tunneling amplitude to 
compensate the S/T energy gap A, so that the Kondo 
effect manifests itself as a zero-bias anomaly (ZBA). An- 
other mechanism of such compensation was discussed in 
Ref. 01, where the energy deficit was covered by the 
conduction electron acceleration at finite bias. In that 
case, the Kondo regime arises in non-equilibrium con- 
dition, so it is fragile against dephasing effects [lil H^. 
In our case the system remains in thermodynamic equi- 
librium around ZBA, so that the only limiting factor of 
this sort is the lifetime of local vibration mode, which is 
usually long enough in comparison with h/Tx- 

To conclude, we demonstrated in this work that 
phonon emission/absorption can induce Kondo tunnel- 
ing in a transition-metal organic complex with even elec- 
tron occupation and a spin singlet ground state, when 
the conventional Kondo effect is suppressed. Unlike the 
situation studied in the current literature 0, 0, 
where the influence of real phonon excitations on the con- 
ventional Kondo effect is discussed, and various kinds of 
side-band satellites due to the polaronic effect are con- 
sidered, we appeal to virtual phonon excitation, so that 
the system remains in a quasi-elastic tunneling regime. 
One of essential ingredients of our theory is that we use 
the dynamical symmetry of the TMOC, which charac- 
terizes both the spin algebra of localized spin itself and 
transitions between various levels of different spin multi- 
plets Il7]. In our case the dynamical symmetry group is 
S0{4. 

Since the tunnel contact between magnetic ion and 
metallic reservoir in TMOC is mediated by a ligand cage, 
one may be sure that the relevant vibration excitations 
are the local phonons characterizing this cage. Although 
we confined ourselves to a specific model with two elec- 
trons in an e subshell, the mechanism is quite general. 
The theory may be easily modified for any system with 
the same structure of the lowest spin multiplets. One 
of such examples is the so called "Fulde molecule" 
schematically representing the spectra of lanthanocenes, 
where the rare-earth magnetic ion is sandwiched between 
two rings of CH radicals. In that case the number of 
electrons in the cage and in the 4/-shell is odd. Another 



candidate is the endofullerene family with atoms [23 or 
magnetic ions within a carbon cage |24l |. The theory of 
vibration-induced Kondo effect may be also generalized 
for the case of degenerate modes and for more compli- 
cated spin- multiplets including half-integer spins. These 
issues will be discussed in forthcoming publications. 
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